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Abstract: The fet and Cambridge/ Aachen inclusive jet finding algorithms for hadron-hadron col- 
lisions can be seen as belonging to a broader class of sequential recombination jet algorithms, 
parametrised by the power of the energy scale in the distance measure. We examine some prop- 
erties of a new member of this class, for which the power is negative. This "anti-fct" algorithm 
essentially behaves like an idealised cone algorithm, in that jets with only soft fragmentation are 
conical, active and passive areas are equal, the area anomalous dimensions are zero, the non-global 
4h ■ logarithms are those of a rigid boundary and the Milan factor is universal. None of these properties 
hold for existing sequential recombination algorithms, nor for cone algorithms with split-merge 
steps, such as SISCone. They are however the identifying characteristics of the collinear unsafe 
plain "iterative cone" algorithm, for which the anti-A^ algorithm provides a natural, fast, infrared 
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Jet clustering algorithms are among the main tools for analysing data from hadronic collisions. 
Their widespread use at the Tevatron and the prospect of unprecedented final-state complexity 
at the upcoming Large Hadron Collider (LHC) have stimulated considerable debate concerning 
the merits of different kinds of jet algorithm. Part of the discussion has centred on the relative 
advantages of sequential recombination (kt p] and Cambridge/ Aachen [2]) and cone (e.g. [3]) jet 
algorithms, with an issue of particular interest being that of the regularity of the boundaries of the 
resulting jets. This is related to the question of their sensitivity to non-perturbative effects like 
hadronisation and underlying event contamination and arises also in the context of experimental 
calibration. 

Recently [I], tools have been developed that allow one, for the first time, to support the debate 
with analytical calculations of the contrasting properties of boundaries of jets within different 
algorithms. One of the main results of that work is that all known infrared and collinear (IRC) 
safe algorithms have the property that soft radiation can provoke irregularities in the boundaries of 
the final jets. This is the case even for SISCone [5], an IRC-safe jet algorithm based on the search 
for stable cones, together with a split-merge step that disentangles overlapping stable cones. One 
might describe current IRC-safe algorithms in general as having a 'soft-adaptable' boundary. 
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A priori it is not clear whether it is better to have regular ('soft-resilient') or less regular (soft- 
adaptable) jets. In particular, regularity implies a certain rigidity in the jet algorithm's ability to 
adapt a jet to the successive branching nature of QCD radiation. On the other hand knowledge 
of the typical shape of jets is often quoted as facilitating experimental calibration of jets, and 
soft-resilience can simplify certain theoretical calculations, as well as eliminate some parts of the 
momentum-resolution loss caused by underlying-event and pileup contamination. 

Examples of jet algorithms with a soft-resilient boundary are the plain "iterative cone" algo- 
rithm, as used for example in the CMS collaboration [6], and fixed-cone algorithms such as Pythia's 
[7j CellJet. The CMS iterative cone takes the hardest object (particle, calorimeter tower) in the 
event, uses it to seed an iterative process of looking for a stable cone, which is then called a jet. 
It then removes all the particles contained in that jet from the event and repeats the procedure 
with the hardest available remaining seed, again and again until no seeds remain. The fixed-cone 
algorithms are similar, but simply define a jet as the cone around the hardest seed, skipping the 
iterative search for a stable cone. Though simple experimentally, both kinds of algorithm have the 
crucial drawback that if applied at particle level they are collinear unsafe, since the hardest particle 
is easily changed by a quasi-collinear splitting, leading to divergences in higher-order perturbative 
calculations 

In this paper it is not our intention to advocate one or other type of algorithm in the debate 
concerning soft-resilient versus soft-adaptable algorithms. Rather, we feel that this debate can be 
more fruitfully served by proposing a simple, IRC safe, soft-resilient jet algorithm, one that leads 
to jets whose shape is not influenced by soft radiation. To do so, we take a quite non-obvious route, 
because instead of making use of the concept of a stable cone, we start by generalising the existing 
sequential recombination algorithms, kt p] and Cambridge/Aachen [2]. 

As usual, one introduces distances dij between entities (particles, pseudojets) i and j and diB 
between entity i and the beam (B). The (inclusive) clustering proceeds by identifying the smallest 
of the distances and if it is a dij recombining entities i and j, while if it is diB calling i a jet and 
removing it from the list of entities. The distances are recalculated and the procedure repeated 
until no entities are left. 

The extension relative to the kt and Cambridge/ Aachen algorithms lies in our definition of the 
distance measures: 

A?. 

d tj = mm(kff, k 2 t J) , (la) 
diB = kg, (lb) 

where A?- = (y, — yj) 2 + (4>i — 4>j) 2 and kti, y% and (pi are respectively the transverse momentum, 
rapidity and azimuth of particle i. In addition to the usual radius parameter R, we have added a 
parameter p to govern the relative power of the energy versus geometrical (A«) scales. 

For p = 1 one recovers the inclusive kt algorithm. It can be shown in general that for p > 
the behaviour of the jet algorithm with respect to soft radiation is rather similar to that observed 
for the kt algorithm, because what matters is the ordering between particles and for finite A this 
is maintained for all positive values of p. The case of p = is special and it corresponds to the 
inclusive Cambridge/ Aachen algorithm. 

1 This is discussed in the appendix in detail for the iterative cone, and there we also introduce the terminology 
iterative cone with split-merge steps (IC-SM) and iterative cone with progressive removal (IC-PR), so as to distinguish 
the two broad classes of iterative cone algorithms. 
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Negative values of p might at first sight seem pathological. We shall see that they are noti The 
behaviour with respect to soft radiation will be similar for all p < 0, so here we will concentrate on 
p = — 1, and refer to it as the "anti-fcf" jet-clustering algorithm. 

2 Characteristics and properties 
2.1 General behaviour 

The functionality of the anti-fct algorithm can be understood by considering an event with a few well- 
separated hard particles with transverse momenta kti, kt2, • • .and many soft particles. The dn = 
min(l/fc| 1 , l/kf^Af.jR 2 between a hard particle 1 and a soft particle i is exclusively determined by 
the transverse momentum of the hard particle and the An separation. The dij between similarly 
separated soft particles will instead be much larger. Therefore soft particles will tend to cluster with 
hard ones long before they cluster among themselves. If a hard particle has no hard neighbours 
within a distance 2R, then it will simply accumulate all the soft particles within a circle of radius 
R, resulting in a perfectly conical jet. 

If another hard particle 2 is present such that R < A 12 < 2R then there will be two hard jets. It 
is not possible for both to be perfectly conical. If kti 2> k%i then jet 1 will be conical and jet 2 will 
be partly conical, since it will miss the part overlapping with jet 1. Instead if kti = kt2 neither jet 
will be conical and the overlapping part will simply be divided by a straight line equally between 
the two. For a general situation, kti ~ ktz, both cones will be clipped, with the boundary b between 
them defined by ARib/kti = A2b/kt2- 

Similarly one can work out what happens with A12 < R. Here particles 1 and 2 will cluster to 
form a single jet. If kti S> kt2 then it will be a conical jet centred on k%. For kti ~ ^t2 the shape 
will instead be more complex, being the union of cones (radius < R) around each hard particle plus 
a cone (of radius R) centred on the final jet. 

The key feature above is that the soft particles do not modify the shape of the jet, while hard 
particles do. I.e. the jet boundary in this algorithm is resilient with respect to soft radiation, but 
flexible with respect to hard radiationH 

The behaviours of different jet algorithms are illustrated in fig. [H We have taken a parton-level 
event together with ~ 10 4 random soft 'ghost' particles (as in [3]) and then clustered them with 
4 different jet algorithms. For each of the partonic jets, we have shown the region within which 
the random ghosts are clustered into that jet. For the kt and Cambridge/ Aachen algorithms, that 
region depends somewhat on the specific set of ghosts and the jagged borders of the jets are a 
consequence of the randomness of the ghosts — the jet algorithm is adaptive in its response to soft 
particles, and that adaptiveness applies also to the ghosts which take part in the clustering. For 
SIS Cone one sees that single-particle jets are regular (though with a radius R/2 — as pointed out 
in [1]), while composite jets have more varied shapes. Finally with the anti-kt algorithm, the hard 
jets are all circular with a radius R, and only the softer jets have more complex shapes. The pair 
of jets near (j> = 5 and y = 2 provides an interesting example in this respect. The left-hand one 
is much softer than the right-hand one. SISCone (and Cam/Aachen) place the boundary between 

2 Note that, for p < 0, min(fc^f, k^) differs from another possible extension, [min(fefj, fc?j)] p , which can lead to 
strange behaviours. 

3 For comparison, 1C-PR algorithms behave as follows: with R < A12 < 2R, the harder of the two jets will be 
fully conical, while the softer will be clipped regardless of whether p t i and pt2 are similar or disparate scales; with 
A12 < R the jet will be just a circle centred on the final jet. 
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Figure 1: A sample parton-level event (generated with Herwig [8]), together with many random soft 
"ghosts", clustered with four different jets algorithms, illustrating the "active" catchment areas of 
the resulting hard jets. For kt and Cam/Aachen the detailed shapes are in part determined by the 
specific set of ghosts used, and change when the ghosts are modified. 

the jets roughly midway between them. Anti-fct instead generates a circular hard jet, which clips a 
lens-shaped region out of the soft one, leaving behind a crescent. 

The above properties of the anti-fc^ algorithm translate into concrete results for various quanti- 
tative properties of jets, as we outline below. 

2.2 Area-related properties 

The most concrete context in which to quantitatively discuss the properties of jet boundaries for 
different algorithms is in the calculation of jet areas. 

Two definitions were given for jet areas in [Ij: the passive area (a) which measures a jet's 
susceptibility to point-like radiation, and the active area (^4) which measures its susceptibility to 
diffuse radiation. The simplest place to observe the impact of soft resilience is in the passive area for 
a jet consisting of a hard particle pi and a soft one p2, separated by a y — (ft distance A12. In usual 
IRC safe jet algorithms (JA), the passive area ajA,R(Ai2) is irR 2 when A12 = 0, but changes when 
A12 is increased. In contrast, since the boundaries of anti-fci jets are unaffected by soft radiation, 
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o(lPJ) A(1PJ) 


cr(lPJ) E(1PJ) 


d D 


s S 


h 


1 0.81 


0.28 


0.56 0.52 


0.45 0.41 


Cam/ Aachen 


1 0.81 


0.26 


0.08 0.08 


0.24 0.19 


SISCone 


1 1/4 





-0.06 0.12 


0.09 0.07 


anti-fc t 


1 1 












Table 1: A summary of main area results for the anti-fc^ algorithm compared to those for other IRC 
safe algorithms in [3]: the passive (a) and active (A) areas for 1-particle jets (IP J), the magnitude of 
the passive/active area fluctuations (a, £), followed by the coefficients of the respective anomalous 
dimensions {d, D; s, S), in the presence of perturbative QCD radiation. All results are normalised 
to irR 2 , and rounded to two decimal figures. For algorithms other than anti-A^, active-area results 
hold only in the small-i? limit, though finite- R corrections are small. 



their passive area is always independent of A12: 

Oanti-fc t ,fl(Al2) = kR 2 . (2) 

Furthermore, since soft particles only cluster among themselves after all clusterings with hard 
particles have occurred, passive and active areas are identical, another feature that is unique to the 
anti-fct algorithm. Thus, specifically for our energy-ordered two-particle configuration, the active 
area is 

Ainti-fci,.R(Al2) = Oanti-fc t ,fl(^12) = • (3) 

In [4] the fact that ajA,i?(0) 7^ ajA,/?(Ai2) (and similarly for the active area) meant that the jet 
areas acquired an anomalous dimension in their dependence on the jet transverse momentum ptj, 
related to emission of perturbative soft particles: 

( ajAjR ) = vri? 2 + djA,R In -?$°L , (4) 

7t6 a s {Rptj) 

where C\ is the colour factor of the hard particle in the jet, Qq is the non-perturbative scale and 
the coefficient of the anomalous dimension, djA,R, is given by 

r2R j a 

d JA ,R= / -i-(a JA , R (0)-7rR 2 ). (5) 



J0 v 

Obviously in this case <i an ti-fct,.R = and similarly for the active area anomalous dimension coefficient 
Aa,nti-fct,.R- One corollary of this is that anti-fct jet areas can be calculated perturbatively, order by 
order, since they are infrared safe. In this respect we recall that they do deviate from ttR 2 for 
configurations with several neighbouring hard particles. 

A summary of these and other properties of the jet area is given in table [H together with a 
comparison to other IRC safe algorithms as studied in [4 J . Fig. [2] illustrates the distribution of 
areas in dijet events at the LHC (generated with Pythia 6.4, with a transverse momentum cut of 
1 TeV on the 2^2 scattering, considering those among the two hardest jets in each event that 
have \y\ < 2), compared with those for other algorithms, and one sees a near 5-function at an area 
of ttR 2 , to be compared to the broader distributions of other algorithms. Fig. [3] shows the mean jet 
area, together with a band representing the event-by-event fluctuations of the area, as a function of 
the jet pt in gg — > gg events, now generated with Herwig. This highlights the independence of the 
area on the jet pt and, once again, the very small fluctuations in the jet area. In this respect, we 
recall that the impact of the UE and pileup on the momentum resolution for jets is related both to 
the typical value of the area (smaller than kt, similar to Cambridge/ Aachen, larger than SISCone) 
and to the extent of the fluctuations of the area, which are close to zero only for anti-fct . 
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Figure 2: Distribution of areas in dijet events at the LHC. We have generated events with Pythia 6.4 
with &pt,min of 1 TeV. Only the two hardest jets have been kept with a further requirement \y\ < 2. 
The area distribution obtained from anti-kt (scaled by 0.1) is compared to the other algorithms, 
(a) passive area at parton level; (b) active area at hadron level including the underlying event and 
pileup corresponding to high luminosity LHC running. 
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Figure 3: Average jet area in dijet events at the LHC. We have generated events with Herwig 6.5 
and only the two hardest jets with \y\ < 2 have been kept. The curves correspond to the average 
jet area at a given pt. The yellow band around the anti-kt line corresponds to the area fluctuations. 
For clarity, the latter are not shown for the other algorithms. The encapsulated graph is a zoom 
on the anti-kt results. Note that the horizontal scale is lnln(pt/ Aqcd) with Aqcd = 200 MeV. 
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Figure 4: Distribution of back-reaction for the anti-fcj algorithm as compared to k t , Cam- 
bridge/Aachen and SISCone. It is calculated for dijet events simulated with Pythia 6.4 in which the 
two hardest jets have p t > 200 GeV and are both situated at \y\ < 2. The back reaction corresponds 
to the net transverse momentum change of each of the two hardest jets due to the reassignment of 
non-pileup particles when one adds high-luminosity LHC pileup to the event (~ 25 pp interactions 
per bunch crossing). 

Back reaction. Suppose one has a hard scattering event that leads to a set of jets {Ji}- If one 
adds a soft event (UE, pileup) to it and reruns the the algorithm, one will obtain a set of jets 
{Ji} that differ in two respects: soft energy will have been added to each jet, and additionally 
the way the particles from the hard event are distributed into jets may also have changed: even 
if one finds the Ji that is closest to the original J\, the two jets will not contain exactly the same 
subset of particles from the original hard event. This was called "back-reaction" in [3]. If the soft 
particles that are added have a density p of transverse momentum per unit area, then for usual 
sequential recombination algorithms the probability, dP/d\nAp\ , of having a back reaction of 
Ap { t B) is O (a sP /Apl B) ) for Ap[ B) > p@ 

For the anti-fcj algorithm, one can show that the probability of back-reaction is suppressed not 
by the amount of back reaction itself, Ap\ , but by the jet transverse momentum ptj, leading to a 
much smaller effect. The impact of this is illustrated in fig. which shows the back reaction that 
occurs for the hard jet events when adding high luminosity LHC pileup to the event. One sees that 
it is strongly suppressed for the anti-Zct algorithm relative to the others, a characteristic that can 
help reduce the smearing of jets' momenta due to the UE and pileup. 

2.3 Other properties 

Non-global logarithms. Resummations of observables involving boundaries in phase space are 
known to involve 'non-global' logarithms. Examples are the jet-mass distribution [9] and energy flow 
distributions in restricted regions \10\ [TT1 [T2] . In both cases there are all-order single logarithmic 
terms (for a hard scale Q, a™ ln n Q/m for jet-mass distributions, a™ ln n Q/E for energy flow), that 
are related to the impact of emissions outside the boundary that radiate a gluon across the boundary 
and so affect the observable defined only on particles within the boundary. 

4 For smaller values of Af>j S ' the full answer requires a resummation that has yet to be carried out, and in the case 
of SISCone, it depends on the nature of the fluctuations of soft event. 
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It was pointed out in [13] that the structure of these non-global terms is more complex with 
clustering jet algorithms, because the boundary of the jet itself depends on the soft radiation^ 
Furthermore this affects even the contribution from independent emissions [14| which remained 
simple in the case of a rigid boundary. So far these effects have been calculated only for the kt 
algorithm, while for the cone and Cambridge/ Aachen algorithms little is known about the non- 
global logarithms other than that they too involve a subtle interaction between the clustering and 
the non-global resummation. 

Because soft radiation does not affect the boundary of anti-fc^ jets, it is straightforward to see 
that their single- logarithmic non-global terms are identically those associated with ideal cones, 
considerably simplifying their determination. 

Milan factor. A crucial ingredient in analytical studies of non-perturbative effects in event and 
jet-shapes is the Milan factor [15 ] 1161 [T7] . which is the correction that relates calculations made in 
a single soft (almost non-perturbative) gluon approximation to calculations in which the soft gluon 
splits at large angle. A remarkable characteristic of the Milan factor ('universality') is that it is 
identical for all event shapes. This is essentially because for all event shapes, if one takes an event 
with hard momenta pi and soft momenta ki , then the event shape can be approximated as 

V({ Pi }, {h}) = V{{ Pi }) + MOi, {Pi})hi (6) 

{ki} 

where fy is a function specific to the event shape observable V. The key feature is the linearity 
of the second term on the right-hand side of eq. ([6]) (see also [IE])- If ; however, the event shape is 
defined for just a jet (or is simply a characteristic of the jet such as its transverse momentum), then 
one loses the linear dependence on soft momenta, since the question of whether one soft particle 
contributes to the observable is affected by its potential clustering with other soft particles. 

In the case of anti-fct jets, the independence of the jet boundary on the soft particle configuration 
means that the approximation eq. ([6]) does hold and the Milan factor retains its 'universal' value 
(M = 1.49 for 3 active non-perturbative flavours [I~5l fTB] ). 

Speed. A relevant issue in order for a jet algorithm to be useful in practice is the computing time 
required to carry out the clustering. The full class of generalised kt algorithms is amenable to fast 
implementation using the techniques of |19j . with the proviso that for p < 0, the specific manner in 
which particles are clustered triggers a worst-case scenario for the Voronoi-diagram based dynamic 
nearest-neighbour graph determination [20] . This means that asymptotically the algorithm takes 
a time O (iV 3 / 2 ) to cluster N particles (rather than iVlniV for the kt algorithm). However for 
N < 20000 the FastJet implementation [21] in any case uses other strategies, which are insensitive 
to this issue, and the &nti-k t clustering is then as fast as k t clustering. 

2.4 Example application: top reconstruction 

One may wonder whether the unusual soft-resilience of the anti-fct algorithm leads to poorer results 
in phenomenological applications. We have investigated various examples and found that in general 
this is not the case. In figure [5] we illustrate this for top mass reconstruction in LHC tt events, as 
simulated with Pythia [7] , where both the t and the i decay hadronically, according to t — > bW + — > 

5 There are cases where, despite the added analytical complexity, this is an advantage since it can cancel a significant 
part of the non-global logarithms. 
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Figure 5: Top mass reconstruction in Pythia-simulated LHC ti events. Both the t and the t decay 
hadronically, t — > bW + — > bqq and i — > bW~ — * bqq. All jet algorithms have been used with 
R = 0.4. 

bqq and i — > rjVF - — » The following simple analysis procedure has been used: we select events 
with at least 6 hard jets (with p t above 10 GeV and \y\ < 5); we assume that both 6-jets have been 
tagged; the 4 hardest remaining jets are paired according to the combination that better reproduces 
the W masses; finally, the W- and 6-jets are recombined to minimise the mass-difference between the 
two i-jets. We use the same four algorithms shown in figure [H now with R = 0.4, and find that they 
all behave rather similarly, with Cambridge/ Aachen and anti-A^ performing marginally better than 
the other two. We note that the difference between various choices of R can be substantially greater 
than the differences between the various algorithms at a given R. One should also bear in mind that 
top reconstruction, near threshold, with a moderate jet radius and no pileup is a relatively simple 
application for inclusive jet algorithms. One expects that greater differences between algorithms 
will be seen in other contexts, an extreme example being multi-scale processes with large pileup. 

3 Conclusions 

There starts to be a certain choice of infrared and collinear safe inclusive jet algorithms for hadron 
colliders. As we have seen, some of these (kt and Cambridge/ Aachen) belong to a more general 
class of sequential recombination algorithms, parametrised by a continuous parameter p, which sets 
the power of the transverse momentum scale relative to the geometrical distance (p = 1 gives kt, 
p = gives Cambridge/ Aachen). 

Rather surprisingly, taking p to be negative also yields an algorithm that is infrared and collinear 
safe and has sensible phenomenological behaviour. We have specifically studied p = —1, the 
"anti-fct" algorithm, and highlighted various simple theoretical properties, notably the resilience 
of its jet boundaries with respect to soft radiation. The other properties that we've discussed are 
essentially consequences of this feature. These properties are characteristic also of certain "iterative 
cone" (and fixed-cone) algorithms, those with progressive removal (IC-PR) of the stable cones, as 
used for example in CMS. However in the anti-^ algorithm these properties are obtained without 
having to pay the price of collinear unsafety. Therefore the anti-fc^ algorithm should be a good 
candidate as a replacement algorithm for IC-PR algorithms. 

One might worry that the resilience (or rigidity) with respect to soft radiation could worsen 



Pythia6.4 SISCone (f=0.75) 




LHC (no pileup) top mass = 1 75 GeV 



9 



the practical performance of the tool. This seems not to be the case, at the least in examples 
examined so far, including the one shown in section 12.41 top reconstruction. One reason for this 
might be that any loss of resolution due to inflexibility in adapting to perturbative branching may 
be compensated for by the reduced fluctuations caused by the underlying event radiation (due 
to suppression of area fluctuations and back reaction). It should also be kept in mind that soft 
particles contribute only a modest component of the overall jet momentum, and the algorithm 
remains flexible in its adaptation to hard (sub)structure in the jets. In this respect it might also 
be interesting in future phenomenological investigations to examine less negative values of p, for 
which the 'hard'-adaptability extends down to softer momenta than is the case with p = —1. 



Note. Subsequent to the first presentation of the anti-fcj algorithm and its main properties at the 
June 2007 Les Houches workshop on Physics at TeV colliders, it was brought to our attention by 
Pierre- Antoine Delsart and Peter Loch that a variant of the kt sequential recombination algorithm 
had been coded within the ATLAS software framework and called reverse-^. It has distance 
measures d{j = max(fc^, A^-)i? 2 /A?- and d{B = k^, and recombines the pair of objects with the 
largest dy unless a diB is larger in which case % is called a jet. We observe that these distance 
measures are just the reciprocals of those for anti-fct; together with the recombination of the most 
distant pair first, this causes the algorithm to produce identical jets to anti-fcf. 
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Appendix: comment on the "iterative cone" 

There are two broad classes of iterative cone algorithms: some find stables cones (of radius R) by 
iterating from all seeds (and sometimes midpoint seeds) and then run some split-merge procedure 
to deal with overlapping stable cones. These algorithms, iterative cones with a split-merge step 
(IC-SM), are used at the Tevatron (for instance the MidPoint algorithm) and it is our understanding 
that the ATLAS cone jet algorithm is of this type too. A second class takes the hardest seed particle 
in the event, iterates to find a stable cone, calls it a jet, removes all particles contained in that jet 
from the event, and then repeats the procedure with the remaining particles, over and over until 
there remain no seeds. This iterative-cone with progressive removal of particles (IC-PR) is the 
iterative cone in the CMS experiment]^ 

IC-SM algorithms have been widely studied, are known to suffer from infrared (IR) safety issues 
(see e.g. [3l [231 E] ) - The IR safe equivalent is SISCone. Though referred to as "cones," they do not 
as a rule give jets with an area of nR 2 [3]. 

6 It is often ascribed to UA1, ref. [22], however the algorithm described there is not based on the iteration of cones. 
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IC-PR algorithms have received less attention, and they behave rather differently from the IC- 
SM variety. For reference we therefore document some aspects of them here. A main point to note 
is that they suffer from a collinear safety issue (in the limit of a fine calorimeter). This can be 
illustrated with the following configuration for an algorithm with radius R: 

Pl : y = 0, <j) = 0.0R, Pt = 130 GeV , (7a) 
p 2 : y = 0, <f, = 0.7R, p t = 200 GeV , (7b) 
p 3 : y = 0, <p= 1.5R, p t = 90 GeV. (7c) 

The hardest seed is P2, and in the small R limit (in which the results are independent of the 
recombination scheme) it leads to a jet that contains all particles, and is centred at eft ~ 0.65-R. If 
particle p2 splits collinearly into 

P2a : V = 0, = 0.7R, p t = 120 GeV , (8a) 
P2b ■ V = 0, (f> = 0.7R, p t = 80 GeV , (8b) 

then pi is now the hardest seed and it leads to a jet at (j) ~ 0A2R which contains p±, p2 a , P2b- That 
leaves p%, which seeds a second jet, and so the number of hard jets is modified by the collinear 
splitting. This will for example lead to divergences at NNNLO in inclusive jet cross sections, at 
NNLO in 3-jet and W+2-jet cross sections and at NLO in jet-mass distributions for 3-jet and W+2- 
jet events. We recall, as discussed e.g. in [5], that a cross section that diverges at N p LO can only 
be reliably calculated up to N p_2 LO. In this respect IC-PR are 'dangerous' to the same extent as 
midpoint-variants of the IC-SM algorithms^ 

Unlike their IC-SM cousins, and the collinear-safety issue aside, IC-PR algorithms do have all 
the features of an "ideal" cone algorithm as described here for the anti-fc^ algorithm, in particular 
they are soft-resilient and give circular jets of radius R. Thus the anti-A^ algorithm seems a natural 
replacement for them, especially as they happen to be identical perturbatively at NLO in the 
inclusive jet spectrum (if used with the same recombination scheme). An alternative replacement 
would be the following: find all stable cones (as in SISCone), label the hardest one a jet, remove its 
particles from the event, and then repeat the procedure on the remaining particles, over and over. 
It too would have the property of soft-resilience, but would differ perturbatively from IC-PR at 
NLO. This, together with the issue of computational speed, leads us to prefer the anti-fc^ algorithm 
as a replacement. 
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